NASA Technical Memorandum 83753 


N84-30273 


Orbital Stability in Combined Uniform Axial 
and Three-Dimensional Wiggler Magnetic 
Fields for Free-Electron Lasers 


Shayne Johnston 
Lewis Research Center 
Cleveland, Ohio 


August 1984 





ORBITAL STABILITY IN COMBINED UNIFORM AXIAL AND THREE-DIMENSIONAL 
WIGGLER MAGNETIC FIELDS FOR FREE-ELECTRON LASERS 

R. Shayne Johnston* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

Zachary [Phys. Rev. A 2_9 (6)m 3224 (1984)] has recently analyzed 
the instability of relativistic-electron helical trajectories in com- 
bined uniform axial and helical wiggler magnetic fields when the 
radial variation of the wiggler field is taken into account. It is 
shown here that the type II instability comprised of secular terms 
growing linearly in time, identified by Zachary and earlier by 
Diament [Phys. Rev. A 23 (5), 2537 (1981)], is an artifact of simple 
perturbation theory. A multiple-time-scale perturbation analysis 
reveals a nonsecular evolution on a slower time scale which accommo- 
dates an arbitrary initial perturbation. It is shown that, in the 
absence of exponential instability, the electron seeks a modified 
helical orbit more appropriate to its perturbed state and oscillates 
stably about it. Thus, the perturbed motion is oscillatory but 
nonsecular, and hence the helical orbits are stable. 
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I. INTRODUCTION 

A recent paper by W. W. Zachary^" extends to arbitrary radius 

2 

some earlier work by P. Diament concerning the instability of rela- 
tivistic-electron helical trajectories in combined uniform axial and 
helical wiggler magnetic fields when the radial variation and axial 
component of the wiggler field are taken into account. These authors 
find in particular that, in the absence of exponential instability, 
there remains a weak instability comprised of secular terms which grow 
linearly in time, casting doubt on the suitability of these helical 
orbits as a basis for free-electron laser calculations. Here, it is 
shown that a refined perturbation analysis employing a multiple-time- 
scale formalism reveals a nonsecular evolution of the orbit on a 
slower time scale which accommodates arbitrary perturbing initial 
conditions. In the absence of exponential instability, the electron 
seeks a modified helical orbit more appropriate to its perturbed state 
and oscillates stably about this modified orbit. Thus, the secular 
terms can be said to manifest physical stability of the orbits rather 
than instability. 

The plan of the paper is as follows. In £ec. II, the formalism 
of the multiple-time-scale perturbation theory is introduced. In 
Sec, iil, the first-order nonsecular solution is obtained by Laplace 
transformation methods and it is shown how the parameters of the zero- 
order helical orbit are modified by the perturbing initial conditions. 
In Sec. IV, the analogous parameter shifts in response to a small 
change in the guide magnetic field. are derived. Finally in Sec. V, 
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2 

the weak instability of Diament is interpreted in this context. 


II. MULTIPLE-TIME-SCALE PERTURBATION THEORY 


Zachary"^ considers the trajectory of a relativistic electron in 
the magnetic-field combination 


^ 2l^ (kr) 

B(r,0,zj = v [ b q z + B w r — — cos (kz - 0)J 


(2-1) 


where B and B denote the magnitudes of the axial guide and helical 
o w 

Wiggler magnetic fields, respectively, k is the wave number of the 
axially periodic wiggler field, 1^ is the modified Bessel function of 
the first kind of order one, and (r,$rf,z) denote cylindrical coordinates 
in space. We adopt Zachary's notation in the following but refine his 
perturbation analysis. 

The equations governing the motion of a relativistic electron 
in the magnetostatic field (2-1) can be written schematically in the 
form 


-* . 
x = y (x , . 


" x 6 } 


(2-2) 


where x denotes the six-dimensional phase-space vector 


-> 

x 


-I r < *■ 2 ' V r' V 


(2-3) 


and (V r ,V ,V^) 316 the cylindrical components of the electron velocity 
multiplied by , the usual relativistic factor. Expressions for the 
six components of are given in Ref. 1. The equations (2-2) admit 
the steady-state helical-orbit solution 
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x* o (t) = £ a, H ^"kut + © i TT/2, K ^ut + k ^©, 0, kau, uJ T , (2-4) 

where © is an arbitrary phase and the parameters a, u and are related 
by the conditions 

y-iku =JC1 ± 2jn w I i (ka) t 1 + < ka) " 2 ] - (2-5) 

o 

\ 

* _1 ku = kc{ (1 - y" 2 ) [ 1 + (ka) 2 ] _1 } 35 , (2-6) 

1 

with jQ. = j e j b / (mctf). 
o,w ' 1 o,w 7 

' 

Consider now perturbations w(t) of the steady-state helical orbit 
X Q (t) by writing 

x(t) = x (t) + £ w(t) + 0(6 2 ). (2-7) 

o 

In other words, we seek a perturbation solution of Eq. (2-2) 
correct to first order in €. where € denotes the order of smallness of 
the perturbing initial conditions. However, to cope with the secular 
terms which arise in a conventional perturbation analysis^, we intro- 
duce a multiple-time-scale formalism. The essential idea of this method 
is to extend the number of independent time variables to remove 
secularities order by order in the perturbation solution; the removal 
of time secularities on a fast time scale determines the time develop- 
ment of the motion on a slower time scale. 1 The solution x(t) is 
expanded in the small parameter € in the form 

x{t) = x^tt^t^t^ + € wttQft^t^, ...) + ... , (2-8) 

where t ,t, ,t„, ... denote a hierarchy of successively slower time 
o 1 2 

scales; 
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dt dt.. 

_2. = i _A 

dt ' dt 



(2-9) 


Operationally, we treat t Q ,t^,t 2 , *** 33 independent variables, 
expanding the time derivative in Eq. (2-2) according to 


d 

dt 




+ € 


+ e F + 


(2-10) 


^1 

Since the number of time variables has been increased, there is 
additional freedom in the perturbation analysis which we make use of 
to remove, order by order, any time secularities which occur in the 
solution. This approach ensures that the perturbation solution 
represented by ( 2 - 8 ) is uniformly valid, order by order. It should be 
emphasized that there is not sufficient freedom in a conventional 
perturbation analysis of Eq. (2-2) to remove secularities order by 
order; only if the conventional expansion procedure is carried out to 
all orders in £ and the secular terms summed, can the slow evolution 
of the motion be determined. However, in a multiple-time-scale 
analysis of Eq. (2-2) , the condition that the first-order solution 


w(t ,t , ...) be nonsecular as t. 
ox L 


determines the slow evolution 


of x q ( t 0 /t^, . . .) on the t time scale. When the multiple-time-scale 
perturbation solution has been obtained to the desired accuracy, one 


returns to the physical variable t by making the replacements t = t, 


2 -> 
t^ = £t, t 2 = € t, ... in the final expressions for x q (t Q ,t^, . . .) , 

vf(t ,t_ ,...) , etc. 
o 1 

Accordingly, our procedure now is to substitute (2-8) and (2-10) 
into Eq. (2-2) , and then to equate to zero the coefficients of 
successive powers of € . To lowest order ( £ °) , we have 
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it 


y(x 


ol 


/ • • • / 


X o6> 


( 2 - 11 ) 


which admits the axisymmetric helical motion (2-4) for a certain class 

of initial conditions . We assume here that the given initial conditions 

fall into this class with perturbations to these initial conditions 

to be considered in next order. Note that the time variable t in (2-4) 

is now to be replaced by t and, furthermore, that the parameters 

a, u and 0, which characterize the helical solution (2-4), are constants 

on the fast time scale t but can vary on the slower time scale t, . 

o 1 

Proceeding to first order in € , we obtain the equation 


i w 


at 


^i 


A i i? 


( 2 - 12 ) 


where A is a matrix whose entries are constant on the fast time scale 
t and which has the explicit form 


A = 


ku 

Ya 

0 

oC 


- JC 

k 

ICO 


-Ka|C 


,-l 


0 


*■ 

ka 

'V 


(5(a) 


-1 


o 

-1 


- t 


, (2-13) 


where we have followed Zachary 1 in defining the parameters 


* = -kujjl o - 2 12 w [ I 1 (ka) + I 1 , (ka)(ka + (ka) -1 )] J , (2-14) 
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$ = Sl Q - 2X1^ (ka) [ 1 + 2 (ka) -2 ] 


(2-14) 


< = + 2/2. ku I* (ka) 

w 1 

* = - 2H w (ka) -1 I 1 (ka) 


The term ^ x^/ ^ t ^ in (2-12) is absent in a conventional pertur- 
bation analysis. It is to be obtained from the lowest-order solution 
(2-4) by differentiating the parameters a, u and 0 with respect to 
t^, and it has the form 


as 


= P + t o q 


(2-15) 


where the vectors p f q are constants on the fast time scale t Q . 
Noting that the time dependences of uft^) and a(t^) are related by 
the equilibrium constraint (2-5) which implies 


ku ^u 

it. 


( + ku JT) 

ka 


a 

at. 


(2-16) 


wo obtain for p and q the expressions 

= da/dt.^ £ 1, O, O, 6, (ku + katfrj ), 2f^] T 

+ k -1 de/dt 1 [ o, ft, l, o, o, o] T 
<1= ^ da/d^ [o, k, 1, 0, 0, o] T , 

where we have defined = - (kuka) 1 ( o( + ku £ ) . The combination of 
parameters ( o< + ku £ ) appearing in (2-16) - (2-18) has been analyzed 
by Zachary 1 , who notes that it can only vanish if (ka) takes on one and 
only one critical value, viz., ka = 0.850. 


(2-17) 


(2-18) 
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III. FIRST-ORDER NONSECULAR SOLUTION 

Combining (2-12) and (2-15) , we seek a nonsecular solution of 
the first-order equation 


d w . , -> . -> 

rr - * . * - P - t„q 


(3-1) 


Equation (3-1) is conveniently solved by Laplace transformation in the 

“V 

variable t in terms of specified perturbing initial conditions w q 

at t =0. Upon introduction of the Laplace transform of w(t )> 
o o 


0 J(s) 


■r 


-St 


dt e 
o 


° 


(3-2) 


there results from (3-1) the algebraic equation 


(A - si) . Ol)(s) = 


-> ->. 2 
-w q + p/s + q/s 


(3-3) 


where 3^ denotes the 6x6 identity matrix. Solving (3-3) forco(s),we 
find 

R(sO . (- W + p/s + q/s ) 

Oi(s) = , (3-4) 

det (A - si) 

where R(s) denotes the transpose of the matrix of cofactors for 
(A - si) . Explicit expressions for the thirty-six entries in R(s) 
are given in the Appendix. 

To evaluate the Laplace inversion integral corresponding to (3-4) 
by means of the Residue Theorem, we require the zeroes of the 
denominator , i . e . , the eigenvalues of the matrix A. Direct evaluation 
of det (A - si) from (2-13) shows that the six eigenvalues are the 
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six roots of the equation 

det (A - si) = s 2 (s^ + bs 2 + d) = 0 , (3-5) 

where, in agreement with Zachary'*', 

b = t [t + (ka) _1 r] + tf _1 ic[ka + (ka)" 1 ] - S" 1 * , (3-6) 

d = - # " 2 ka K { [l + (ka)" 2 ] ( o< + ku t ) - Y -1 k 2 u 2 ] . (3-7) 

As an independent check on the algebra, we have related the results 

3 4 

(3-6) and (3-7) to those obtained by Freund and co-workers ' without 
Laplace transformation. We have, by judicious use of the equilibrium 
condition (2-5) , carefully verified the relations 

b =Sll + Sl\ , d = A 2 A 2 , (3-8) 

where £1 and n 2 are identical to the (possibly complex) character- 
istic response frequencies of the equilibrium helical orbit derived 

3 

by Freund, Johnston and Spr angle and reported again by Freund and 
, 4 

Ganguly , viz., 

-n-1,2 - 1 <“l + “2> 1 7 + 4 “' - 13-9) 

where 

Oi) 2 = ^" 2 k 2 u 2 + V _1 ku 2jl w I 2 (ka) [ ka + (ka)" 1 ] , (3-10) 

W 2 “[-^o " ^" lku ~ 2jI w I l (ka) ] iSl o " ^" lku) 

- y 1 ku 2 _ft w I 2 (ka) [ ka + (ka) 1 ] 

A = - ca" 1 (SI - 2 Y _1 ku) 
o 


(3-12) 
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B = - 2 il w u (c ^ ) _1 j [n. Q - y -1 ku (1 - k 2 a 2 )'] I 2 (ka) 


+ * _1 ku [ (1 + k 2 a 2 ) IJ (ka) - (ka )" 1 (ka)] ]• . (3-13) 


2 2 


-1 


The' conditions for exponential instability of the orbit are, as 
1 2 

shown by Zachary , either d < 0 with (b - 4cl) > 0 , or d > 0 with 
2 

(b - 4d) < 0. The marginal stability condition d « 0 agrees exactly 

4 

with that found by Freund and Ganguly . In the domain d > 0 with 

2 2 2 

(b - 4d) >0, the quantities and XI ^ are real and positive and so 

there is no exponential instability of the orbit. However, Zachary'*' 
finds a non-exponential instability instead in this case, associated 
with secular terms which grow linearly in time. The proper inter- 
pretation of these secular terms in terms of orbital stability is the 
subject of this paper. 

Insertion of (3-5) and (3-8) into (3-4) yields the relation 
R(s) . (- W Q + p/s + q/s 2 ) 

W (s) = — — = — — 5 5 , (3-14) 

s 2 (s 2 +Jl 2 ) <s 2 +SL) 

1 X 


and it remains simply to invert the Laplace transformation to obtain 

w(t ) . The singularities of ci)(s) occur at the simple poles s = - iXt- , 
O 1 

s = ± m 2 , and also at s = 0. It appears from the form of the term 
in (3-14) that the pole at s B 0 may be as high as fourth order. 
However, explicit evaluation of the q term reveals that R(s) . <3 is 
proportional to a and so the pole at s = 0 is at most of third order. 
Accordingly, we rewrite (3-14) in the form 


40 (s) = 


N (s) 

~ 3 2 2 2 2 

s j (s 2 +/1) (S 2 +JT) 

I 2 


(3-15) 
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where the vector numerator N(s) is analytic at s = 0 (as well as at 
the other zeroes of the denominator) and is given by the expression 


N(s) = 


R(s) . (-sw q + p + q/s) 


(3-16) 


The Laplace inversion integral can now be evaluated using the 

2 

Residue Theorem. In the domain d>0 with (b - 4d)>0, the contributions 

from the simple poles at s = - iX).^ an< ^ s = “ iA^ lead to stable 

sinusoidal oscillations at the characteristic frequencies and 

XI . Thus, the solution \*(t ) can be decomposed into an oscillatory 
z o 

part, \T (t ) , and the contribution from the pole ats = 0, vT. (t ), 
osc o * o 

which includes any time secular it ies : 


w (t ) = w (t ) + w (t ) 
o osc o ’ * o 


(3-17) 


Use of the Residue Theorem leads directly to the oscillatory terms 


w 

osc 


<v 


1 

2 <a 2 -n 2 2 > 



N (i ii!^ 



1 



if (1 m 2 ) 


+ 


e 


iil 2 t 



(3-18) 


It remains, then, to investigate the contribution w (t ) from the pole 
at s = 0. 

Since it is of third order, the pole at s = 0 leads to a contri- 
bution of the form 


w, (t ) 



(3-19) 
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where y[t and £ are vectors, constant on the fast time scale t Q . 

Following the philosophy of the multiple-time-scale approach, we seek 

to eliminate the secular time-growing terms in (3-19), i.e., to arrange 
■ ■ -s* 

that ym = 0 and $ = 0 by an appropriate choice of the slow time 
dependence of a(t^)aj>d e(t^). The term J represents a small net phase- 
space displacement of the orbit which need not vanish. 

Upon use of the Residue Theorem, there result the following 
expressions for the vectors jx, : 


s - 


-Clj 


lim 

s*0 




^2 d_|(s) _ 2 (a 2 +fl 2 } 

ds 


N(s)^ 


,(3-20) 


A 


si 2 si 2 

Jt l U 7, 


lim 

s*0 


difrs) 

dS 


(3-21) 




lim 

s-»0 


N(S) 


(3-22) 


Examining (3-20) - (3-22) , we see from (3-16) that we must determine 
the matrices R (o) , (dR/ds)|s = 0 and (d R/ds )|s = 0. The explicit 
formulas for the entries in the cofactor matrix R(s) given in the 
Appendix make these determinations straightforward. 

Consider first the vector ^given by (3-22). From (3-16) , we 
see that 


lim 

s^O 


dR 


N(s) 


= — (°) p x+ dT L = o • 


(3-23) 


From expressions (2-17) , (2-18) for p, q, there results then the 


formula 


- 12 - 


7 = da 2 ( <X + ku 5 ) r , 

* dt^ (ku) (ka) l 0 ' k ' lf °' °' °] 


Since ( o( + kujT) ^ 0 except for the one critical value of ka, 

therefore set 

da . 

^ ‘ ° 

in order to eliminate the quadratic time secularity in (3-19) . 
follows from the constraint (2-16) that we must also have 


du 

dt. 


= 0 


and from (2-18) that q = 0. 

Proceeding to the evaluation of jx and noting that 


lam dN (s) „ , , -> — 

3*0 ds — o ds 


s = 0 


we find, making use of (3-25) , the explicit result 




= (" * 3 


( « + taj) (ka ”s5 * "oS 1 
# 3 ka 


JTL 2 Si 2 

1 2 


_ i de_ > 

k dt i / 


x [0, k, 1, 0, 0, 0] 


Note that the condition that ym. vanish is again a single scalar 
equation rather than a vector equation, this time leading to a 
determination of the slow precession of the phase 9(t^). From 
the condition for a nonsecular solution with JU = 0 is 


(3-24) 


we must 


(3-25) 

It 


(3-26) 


(3-27) 


(3-28) 


(3-28) , 
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de K ( qC + kuX ) 

dt. v 3 _ 2 p. 2 

1 * a 


<ka "o5 + "o 6 > 


(3-29) 


It remains, finally, to evaluate and interpret the constant 
shift vector g in (3-19) . Using the relation 


2 -> 

lim 1 d N(s) 


1 d ^ 


s+0 2 


ds 


2 , 2 
ds 


dR 


P - 


s=0 


ds 


s=0 


w 

o 


(3-30) 


we obtain for the six components of the results 


C = K (-T- jP/ka) 

1 f sil n 2 % 


(ka w o5 + w o6> 


, -Kfc 

* 2 ‘ 3 ‘ Sl\ SL\ 


f <k a w ( 


ku 

+ i. w 

o2 o4 


kjl> +k " (S + kat)]w o3 j 


s 4 -° 


(3-31) 


x- - 1C (u+ + o(a) ,, 

S r = — V — — (ka w„ c + w^ c ) 




o5 06 


^ - K(e<+ kuJT ) „ , 

£ r = — ~ r 7 (ka w + w ) 

6 * kaJl^A^ 5 


These results can be understood in terms of shifts <Sa, Su, S© in the 
parameters a, u, 9 which characterize the equilibrium helical orbit 
(2-4). They admit the interpretation 


= Sa 


s 2 = k S 3 =Se , S 4 - 0 


$ 5 = (ku) Sa + (ka) Su 


S 6 = ^ 


(3-32) 
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where 


C lUcu 

oa = — - 


Kku 

* 3 


(ka W q5 + W 06 } 

si 3 n? 


Su = 5- 22 - 

dt. 


(3-33) 


s© = 




{ • II I W o3 - £ [*< a<W o2 - * ”o3> + T ” 04 ]} 


with d0/dt^ given by formula (3-29) . It follows that the perturbed 
motion (3-17) comprises stable oscillations at frequencies 
about a modified helical orbit whose parameters are shifted slightly 
by the perturbing initial conditions. 

To conclude this section, it should be noted that the quantity 

(ka w + w ) appearing in (3-33) vanishes for an initial perturbation 

05 06 

which conserves energy [" see the equilibrium solution (2-4)^ . Since 
the basic equations (2-2) conserve energy, such an initial energy 
perturbation 6* remains constant and could be incorporated into the 
zero order energy £ . For energy- conserving initial perturbations 
with 0., the secular terms vanish identically and the only 

non-vanishing shift in (3-33) is £ 9. However, this is not the case 
for the response to perturbations in the magnetic field structure, as 
the next two sections will demonstrate. 
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IV. ORBITAL RESPONSE TO A SMALL CHANGE 
IN THE GUIDE MAGNETIC FIELD 


It has recently been proposed that the efficiency of free-electron 
lasers can be enhanced by an appropriate taper in the axial guide 
magnetic field strength. In such a scheme, it is important to under- 
stand the effects of the changes in the guide field on the motion of 
the electrons. In particular, one can ask whether the helical orbits 
persist or are destroyed (by secular terms) as such changes occur. In 
this section, we show by means of a calculation completly analogous 
to that in Sec. Ill that the effects of a small change in the guide 
majgnetic field are two-fold: the parameters of the helical orbit are 

shifted slightly and oscillations about this modified orbit occur. 

In addition to the reason stated above , these results are also important 
to the extent that they reflect a generic resiliency of the helical 
orbits in response to more general magnetic perturbations. 

Consider a small change £b in the guide magnetic field where 

o 

Kb I « B , and define 
1 o' o 


s_a 

o 


1*1 * B 0 

me V 


(4-1) 


The linearized equations of motion £cf. 


4 

w 


= A . — J? - t q +&n. [ o , 0 , 


(3-1)"] then take the form 
0, -kau, 0, 0] T , (4-2) 


where we treat SB q as a first-order perturbation and where the 
quantities A, p and q are given as before by (2-13), (2-17) and (2-18), 
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respect ive ly . It is clear from (4-2) that the net effect of the 
perturbation £B q on our calculations in Sec. Ill will be to change the 
vector by an amount 



_ i] 

kau SJl [ 0, 0, 0, 1, 0, Oj 


(4-3) 


For simplicity/ we set the initial perturbation w^ = 0 and take (4-3) 
instead as the source of the perturbed motion. 

The solution of (4-2) again takes the form 

v^(t ) = w (t ) + + Mt + "(t 2 /2 , (4-4) 

o osc o * ^ o * o 

where ^gc^o) i s given by (3-18) and jx and ^ by (3-20) - (3-22) , 
provided we now take 


N(s) = R(s) . ( p + %p + q/s) 


(4-5) 


As in Sec. Ill, the condition f = 0 implies da/dt^ = 0 and 
hence du/dt^ = 0 and q = 0. The condition jx - 0 again leads to an 
expression for the slow evolution of the phase 8 (t^) , viz.. 


d8 _ k 2 u 2 k a K rjj 

dt i * 2 tfn 2 n 2 ° 


(4^6) 


Finally, the vector £ can be calculated and interpreted according to 
(3-32), i.e., in terms of shifts Sa, Su and $9 in the parameters 
a, u and 0 which characterize the unperturbed helical orbit. Our 
results for these shifts are 


u 


x 2 ji 2 n 2 


S a 


2 2 

( 1 + k a ) 


(4-7) 
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S Q 
S-flo 


= 0 


( 4 - 8 ) 


£u _ K A (ku ) 2 

^ A 0 y 2 A 2 Jl 2 


( 4 - 9 ) 


where the ratio form displays the rates of change of these quantities 
with respect to a smoothly varying guide field. Our results imply 
that the helical orbits persist in the presence of such a smooth 
variation. 


V. DIAMENT'S WEAK INSTABILITY 
2 

The weak instability of Diament differs somewhat from the type II 
instability of Zachary 1 which was discussed in Sec. III. Diament takes 
for the unperturbed motion the helical orbit associated with the ideal- 
ized wiggler field in which all radial dependence is ignored. He then 
treats the radially-dependent corrections to the wiggler field as a 
perturbation and asks what will be the first-order effect of these 
realizability corrections on the idealized helical motion. As Zachary 

points out 1 , this approach is clearly limited to the case of small ka. 

2 

Diament finds that when the idealized orbit is exponentially stable, 
there remains a weak instability associated with secular terms growing 
linearly in time except for special choices of perturbing initial 
conditions. He concludes from these results that realizable wigglers 
tend to induce outwardly spiralling motion in general. The purpose of 
this section is to properly interpret the secular terms and to show 
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this conclusion to be erroneous. 

2 

Following Diament , we take for the unperturbed magnetic field 
the common approximation 

Bf. = B w cos (kz - 0) , = B w sin (kz - 0) , B^ = B ft , (5—1) 

which is the limit of (2-1) as r -► 0. The radially-dependent corrections 
to (5-1), i.e., the terms which must be added to yield (2-1), are thus 


SB r = B w [2 IJ (kr) - l] 

cos (kz - 

0) 


(5-2) 

= B w[ 2 < kr > _1 1 1 (to) 

- lj sin 

(kz - <f>) 

9 

(5-3) 

SB^ = - 2 B w 1^ (kr) sin 

(kz - 0) 

9 


(5-4) 


and are to be treated as a small perturbation. The limiting field 
(5-1) admits a steady-state helical orbit as a solution of the 
equations of motion (2-2) which again takes the form 

T 

5-5) 

However, relations (2-5) between the parameters a and u now simplifies 
to 

<“ 1 ku =Si - (ka)" 1 -fl-w . (5-6) 

o 

We now proceed with a multiple-time-scale perturbation calculation 
of the first-order correction to the helical motion (5-5) , following 
closely our work .in Secs. Ill and IV. 


t Q (t) = £a , ^"kut + 9 - ^/2 , 5 ( 1 ut + 


k _1 9 


kau 


i 
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The linearized equations of motion take the form [cf. (4-2)] 

— ► *4 —4 . -4 *4 

w = A' . w - p* - -t^q‘ + & , (5-7) 


where henceforth a prime will denote association with the approximation 
(5-1) . The matrix A' has exactly the same structure as the matrix A 
given by (2-13), but the definitions of the parameters , 5 , K , 'P 
must be modified as follows: 

<*' =-* _1 k 2 u 2 , £' = -Sl Q + 2*" 1 ku 

(5-8) 

|C* = + ku , y = - JI W . 

In fact, most of our earlier calculations can be adopted without 
change provided that the replacements (5-8) are made; in particular, 
the formulas given in the Appendix for the entries in the cofactor 
matrix R(s) remain valid. The characteristic response frequencies 
of the equilibrium helical orbit now simplify to 

X " 


JCIJ = tf-^ku , (i^) 2 = ( y _1 ku -Jl 0 )[x" 1 ku -Jl o (1 + k 2 a 2 )] , 

(5-9) 

2 

and the condition for exponential instability reduces to (.Sip i n 
(5-9) being negative. 

The vectors p’ and q 1 in (5-7) again arise from the slow-time-scale 
term ix^/^t^ with (2-17) and (2-18) simplifying to the expressions 


p* ='^da/dt 1 [l, 0, 0, 0,il Q , (-ka)" 1 Mf -1 ku -.fl 0 >] T 
+ k -1 dO/d^ [0, k, 1, 0, 0, 0] T , 


(5-10) 


9 
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q' = - (ka )" 1 -fl Q ) da/d^ [o, k, 1 , 0, 0 , o] T . ( 5 - 11 ) 

The source term & in (5-7) is due to the realizability corrections 
(5-2) - (5-4) and takes the form 


= + uJI w [2 ka (1 + k~ 2 a" 2 ) I (ka) - l] [o, 0, 0, 1, 0, o] T . 

(5-12) 


From the form of (5-7) , we see that the net effect of ^ is again to 
modify the vector p ' , just as in Sec. IV. For simplicity, we again 
set w* o = 0. Since the rest of the calculation proceeds exactly as in 
Sec. IV, we shall simply state and interpret the results. 

The slow phase precession is found to be 


de 

dt. 


ku SL 


- — ^ [2 ka (1 + k“ 2 a~ 2 ) I. (ka) - 1 1 , (5-13) 

1 tfw) 2 L 1 J 


:^) 


with da/dt^ and du/dt^ both zero. The shifts £a, S&, Su in the 
idealized helical orbit parameters (5-5) are as follows: 

Q 2 2 2 

•— — w (1 + k a ) T „ , . ,-2 - 2 V „ „ , ,1 

Sa = 5 - 5 I 2 ka (1 + k a ) I (ka) -1 I , 

(.OI) k a 

(5-14) 


Se = 0 

Su = - 


u ~ [ 2 ka (1 + k" 2 a"? I (ka) -1 1 . 


(5-15) 

(5-16) 


As in the previous sections, the corresponding perturbed motion is now 
nonsecular. 

These results admit the following interpretation when exponential 
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instability is absent. An electron assumed to be executing an 
idealized helical orbit satisfies relation (5-6) , 

a" 1 ku* = SL i (ka') _1 _a 
o w 

whereas, in the presence of the realizability corrections (5-2) 
the proper relation between a and u is (2-5), i.e., 

tf -1 ku = fl - 2SI I. (ka) (1 + k“ 2 a“ 2 ) 
owl 

The electron in question is therefore not on the helical orbit appro- 
priate to its perturbed state and so, in accordance qualitatively with 
our rsults in Secs. Ill and IV, it seeks out this proper helical 
orbit by shifting its orbit parameters and by oscillating stably about 
this modified helix. To support this interpretation, set a = (a 1 + £a) 
and u = (u* + Su) in (5-17) and (5-18) , with &a. and £u given by (5-14) 
and (5-16) , respectively. Upon subtraction of (5-17) from (5-18) , 
the consistency condition is seen to be 

V” 1 k Su - (ka 2 )" 1 !! Sa = - (ka) -1 !! [ 2 ka (1 +k -2 a -2 ) I. (ka) - ll . 

w w l. 1 J 

(5-19) 

To leading order in (ka) , i.e., in the small expansion parameter, the 

dominant contribution to the left-hand side of (5-19) comes from £a and 

condition (5-19) is found to be verified. 

In conclusion, it should be stressed that realizable wigglers do 

not induce outwardly spiralling motion in general. We have shown here 

2 

that the secular terms found by Diament are just an artifact of a 


(5-17) 

- (5-4), 

(5-18) 
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direct perturbation calcillation and that they disappear in a more 

2 

refined analysis. Diament's assertion that, even if the secular 
growth were bounded, the electrons would typically strike the wall of 
the drift tube before turning around, is not correct. The radial 
excursion associated with the shift (5-14) and the amplitude of the 
stable oscillations is of the same order of smallness as the perturba- 
tion parameter (ka) , and the smallness of (ka) is a necessary 
prerequisite for Diament's analysis to be valid. 
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APPENDIX: ENTRIES IN THE COFACTOR MATRIX 

In this Appendix, we summarize the< entries in R(s) , the transpose 
of the matrix of cofactors for (A - si) . Each of the thirty-six 
entries below results from the evaluation of an appropriate five-by-five 
determinant. 

R u = “ s 3 jV + t [*+ (ka )" 1 j] + tf" 1 * [ka + (ka) -1 ] | 

R 12 = s 2 + + (ka)" 1 ^] 


R 13 

- -s 2 ^- 1 ! 

<ka 

[ ¥ + (ka) 

■ l s] • 



R 14 

= - s 2 X _1 

{ 

s 2 + Tf 1 * 

£ ka + (ka) 

*1} • 


R 15 

= - sT 1 

[3 

s 2 - a -1 * 

(t" ka^)] 

• 


R 16 

= SX X [1 

, 2 
’ s 

+ (Vka) _1 

K (H' -ka$ )' 

1 


R 21 

* s 2 (tfa) 

- x [ 

2 

kus + ku 

(+ 2 + ¥" 1 Kfca)+ (ka) _1 t (o( 

+ kuX)J 

R 22 

= -s|s 4 + 

2 

s 

[+(*+ k _1 

a" 1 !) + X" 1 

(kaK -*)] -T 2 

kaK<* J 


R 23 = " sK(ya) -1 [s 2 -V" 1 ^ -y _1 ku (S +kat)] 

R 24 = S < ^ a ) _;L { s 2 [V'^u + (ka)" 1 'V] + ^“ 2 kukaHj 

R 25 = “ (^ a )” 1 (s 2 + V _1 kaK) [ s 2 - V* 1 (tf+kujT)] - { * a) "V + 2 

R 26 “ s2 ( ^ a)-1 [ (ka)” 1 * 2 - y^kut- V” 1 K.*] + (^a)" 1 * ( <* + ku£ ) 

r 31 = s 2 ( ^ _ 1 kuK-^t) 
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R 32 = - sY - 1 K& (s 2 - If ”^c< ) 

R 33 * “ s [ s4 + s 2 k’ 1 a' 1 Jf) + y _1 (k’V’ 1 *-*)] 

- U^ka)" 1 ^ + ku( X + ka'V )] } 

R 34 = - sX 1 (fs 2 - V~ 2 kuK) 

R 35 * ^ it-'K-S* +#~ 3 K(* + ku?) 

R, A * - y -1 S 4 - (ka#) -1 S 2 ( }f" 1 K+^ - 5 f -1 ka o( ) 

36 

+ ( 3 3 ka) “ 1 * ( 0 ( + ku J ) 

R 41 - - s 2 ^s 2 +« _1 o( K t ka '+ (ka)" 1 ]+ Jf^kuK [t+ . 

R42 “ s 3 k" 1 1 C ( ^? + ka t ) 

R43 = - s 3 |C ( T + ka t ) 

R 44 = - s 3 -[ s 2 + V 1 * [ka + (ka)" 1 ] } 

R 45 - -s 2 [ 5 s 2 - * - kaS)] 

R 46 = s2 [^s 2 + ('tfka)" 1 K ( V “ ka^ )] 

R 51 = s3 (ka) 1 < * ♦ “ * ” 1 kuK) 

R 52 = s 2 k~ 1 K<s 2 - X _1 <) 

R 53 = - s 2 K (s 2 - V - 1 o< ) 

R 54 = s 2 (ka)"" 1 < * s 2 - ¥“ 2 kuK) 

R 55 = - s £ s 4 + s 2 ( H* 2 + 3 ^ka.1? - ){ *<*) “V 2 k )C { a + uS* ) J 
R 56 = - s | s 2 [(ka)" 1 ^ 2 + X” 1 *] - ( H^^KWa + uf ) J 
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61 


62 


63 


64 


R 


65 


66 


s 3 < *+ - H '^cuB 

- s 2 ** (s 2 
.Wcs 2 -^* 

- s 2 OVs 2 " ^ _2 kulO 

. s{s n s U'-'' 1,n+r1 "' 

. (jV'iii « ■ + ’“ y ' ) J ‘ 


S>1 • 

• Tf" 1 *] 
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